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Abstract 

We study the properties of the diffuse y-ray background around the Galactic plane at energies 20 - 200 GeV. 
We find that the spectrum of this emission possesses significant spacial variations with respect to the average 
smooth component. The positions and shapes of these spectral features change with the direction on the sky. 
We therefore argue, that the spectral feature around 130 GeV, found in several regions around the Galactic 
Center and in the Galactic plane in HIE El 0L can not be interpreted with confidence as a y-ray line, but may 
be a component of the diffuse background and can be of instrumental or astrophysical origin. Therefore, the 
dark matter origin of this spectral feature becomes dubious. 



1. Introduction 

It has been recently reported in and further investigated in |Ql and in 1 4 j3 that the y-ray emission 
from the region around the Galactic Center (GC) exhibits a line-like excess at the energies ~ 130 GeV. 
An interest to this result is based on the expectation that any signal of astrophysical origin at high energies 
would have a broad (compared to the Fermi spectral resolution) spectral shape. Diffuse emission with the 
line-like spectrum has therefore been considered as an exotic one, e.g. as a "smoking gun" for dark matter 
annihilation (see e.g. for review). In particular, "Higgs in space" scenario [8] predicts a y-ray line 
at ~ 130 GeV for the Higgs mass around 125 GeV as seen by the LHC 

The region of JTt] and H was selected by maximizing signal-to-noise ratio for the expected dark matter 
annihilation signal. The preprints of [1] and [2] were followed by [3J where the claim was confirmed and it 
was demonstrated that a similar excess originates from several regions of the size ~ 3° around the Galactic 
plane. A number of works 03j[lj][l2|[l^[l^[l5j], have discussed possible interpretations and origin of this 



spectral feature, see [ 7 ] for the review. The search for y-ray lines, based on the 2-year data, performed by the 
Fermi collaboration [ 16] did not reveal any lines but had not comment on the origin of the observed excess. 

It was demonstrated in the first version of this paper ( 1205.4700vl) that spectral features with the signifi- 
cance, similar to that of the excess observed in the Galactic center region around 130 GeV can be also found 
at other energies in different regions of the sky. It was then argued in 1 1/7, 18 1 that the only "significant 



feature" is the one in the GC, while all the other features, found in 1205.4700vl can be considered as pure 
statistical fluctuations. 



'The latter paper has appeared after the first version of the present work 1205.4700vl 
2 see also talks by C. Weniger at IDM-2012 OH and COSMO-2012 |2(J] conferences 
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Figure 1: Difference in the number of photons between the energy bin 120-140 GeV and the half-sum of two adjacent bins (100-120 
and 140-160 GeV). The regions with positive and negative excesses around the background are clearly visible. Three most significant 
regions from |3] are shown with white circles. 

We believe that the current data do not allow to reach a definitive conclusion here. Moreover, the most 
relevant question is actually different. To support the DM interpretation of the GC line, one should try to 
determine: 

- Is it possible to exclude the existence of features in several regions around the Galactic plane (apart 
from the GC)? 

We demonstrate in this work that the answer to this question is negative - the data does not allow to separate 
with confidence a feature in the Galactic Center from other features in the Galactic plane. Moreover, common 
interpretation of all the features is possible, providing additional evidence for such a hypothesis. Therefore 
that the nature of all these deviations from simple featureless spectral models at energies 50-200 GeV should 
be discussed, rather then some peculiar properties of the Galactic center region. Until this question is settled 
the DM origin of the 130 GeV line (or a pair of lines at 1 10 and 130 GeV) remains dubious. 



2. Identifying regions with spectral features in the Galactic planes 

For the analysis, presented below we used 209 weeks of Fermi data and v9r23pl Fermi Software. We 
filter the photons with the expression (data_QUAL==1 ) && 

(LAT_CONFIG==l) && (ABS (ROCK_ANGLE) <52 ) && ( ( STOP<3 52 7 7 3 2 . ) II 

(START>352814402 . 0) ) recommended by the Fermi-team in order to exclude a bright solar flare that 
affected Fermi on March 8th, 20120 

2.1. Residual maps 

We start by searching for regions around the Galactic plane (|Z| < 30°, \b\ < 5°) which have excesses at 
different energies. To this end we split the energy range 80-200 GeV into 6 energy bins, 20 GeV each (at 



3 The large area telescope (LAT) on board the Fermi satellite is a pair-conversion gamma-ray detector operating between 20 MeV 
and 300 GeV. The LAT has a wide field of view of ~ 2.4 sr at 1 GeV, and observes the entire sky every two orbits (~3 hr for a Fermi 
orbit at an altitude of ~565 km, full details of the instruments are given in 12111 ). 
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Figure 2: Same as in Fig.[T]but showing the difference of 100-120 GeV and 140-160 GeV energy bins with the half- sum of two 
corresponding adjacent bins. Yellow and magenta circles mark bright excesses on these maps, see text for the details. The radii of all 
circles are 3° 

energies of interest the spectral resolution of Fermi is about 10%). In each energy bin we build a count maps 
using CLEAN class of photons and smooth each map with a tophat filter with the radius of 3°, obtaining 
images I,, i — 1 ... 6. For each image i~, we then define a "background" as 2?,- = h (J;_i + i~,+i). In order to 
identify regions with spectral features we built then a set of images Mi = li/S/ - 1. Such a map for 120- 
140 GeV energy band is shown in Fig.Q] The positive (negative) color coding at these images corresponds 
to the regions with positive (negative) excess in the spectrum. Three regions with the most significant excess 
from J2I] are marked on the map as white circles. The positions of white circles obviously correspond to the 
positions of (some of) the brightest excesses on this map. It is clearly seen in Fig.Qjthat emission around the 
GC region is shifted from (I, b) = (0, 0). 

Notice, however, that excesses of similar significance are visible not only on 120-140 GeV map (Fig.Q]), 
but also on the maps 100-120 GeV and 140-160 GeV (Fig. [2j». We identified three additional regions that 
show high signal/background ratio: Reg 1 and Reg 2 regions are bright on 100-120 GeV map, while Reg 3 
region is bright on 140-160 GeV. All these regions are circles with 3° radius, as are the regions of [3]Q 

The resulting spectra in the energy range 50-200 GeV extracted from the regions Reg 1, REG 2, and 
Reg 3 together with the Central region (coinciding with the region of the most significant excess of J3t] 
around the GC) are shown in Fig. [3] The spectra from all the regions Reg 1, Reg 2, and Reg 3 clearly 
show a number of features at different energies. The thickness of the lines representing the spectra in Fig. [3] 
is defined by the formal lcr statistical error. 
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Table 1: The values of ^ 2 /[d.o.f.] when fitting the power law model (TJ to the data in the energy range 50-200 GeV (column (a)) and 
20-150 GeV energy range (column (b)). Columns (c) and (d) show the chance probability p to get observed x 2 values (as defined 
by Eq. (2))- The significances that correspond to these p-values are given in columns £50-200 ar >d S 20-150- The change of the fitting 
interval may affect the p-values by as much as an order of magnitude. 





Prominent features, 


Observed 


Predicted 


Significance 


Region (1°, b°) 


energy range 


counts 


by model 


of spectral feature 


(a) 


(b) 


(c) 


(d) 


(e) 


Reg 1 (18.72°, 0.57°) 


105-120 GeV 


29 


10.46 


4.63 


Reg 2 (335.76°, 1.255°) 


70-90 GeV 


75 


40.95 


4.72 


Reg 2 (335.76, 1.255° 


95-1 10 GeV 


32 


16.83 


3.23 


Reg 3 (12.22° ,0.20°) 


80-85 GeV 


13 


5.25 


2.74 


Central (359°,-0.7°) 


125-135 GeV 


25 


3.53 


7.34 




105-115 GeV 


17 


5.99 


3.58 



Table 2: The most significant features in the regions REG 1, REG 2, REG 3, and CENTRAL, as well as the features in the region REG 2. 
We fit the background counts to the power law model in the whole 50-200 GeV energy range (unless otherwise specified) and determine 
the most prominent deviations from this simple background model and their significance (column (e)). This significance of the spectral 
feature is defined as the Poisson probability to observe N counts or more, provided that the spectral model predicts A counts (specified 
in the column (rfj)0 



2.2. Analysis of the spectral features 

In order to quantify the significance of the observed spectral features in the regions Reg 1, Reg 2, 
Reg 3 we start with performing a power law fit 

N(E) = N m (E/l00 GeVr r , (1) 

to the data. The reduced^ 2 and the resulting p-values are presented in the Table Q] columns (a-d). Here the 
p-values are defined as 

Pn(x 2 )dx 2 (2) 

4 

where P„(x 2 ) is the probability density function of the^ 2 distribution with n degrees of freedom. 

Looking at the quality of fit to the power law model in 50-200 GeV energy range one could have con- 
cluded that with 20 - 40% probability the data is consistent with being described purely by the power law 
without any features (even in the Central region this probability is more than 3%). Based on a similar 
analysis, it was argued (see e.g. |T^, 2(J) that the features in the regions Reg 1, Reg 2, Reg 3 are not 
statistically significant. 



4 Any difference between the "residual maps" between CLEAN and ULTRACLEAN events would signify that these residuals are 
parts of some known systematics (e.g. contamination from cosmic rays) 0]. We nave verified that one obtains the same maps when 
using ULTRACLEAN photon class. 

"This definition of significance is different from the one, used in (^|. See the discussion in Sec. [3] 
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Table 3: Change in the total x 2 when adding a Gaussian at 130 GeV (CENTRAL), at 112.5 GeV (REG 1) and at 85 GeV (REG 3) as 
compared to the power law fits, shown in Table [TJ A formal significance of these features is defined as yj&x 1 ( as we add only one 
degree of freedom when fitting a normalization of the Gaussian). The features in the region REG 2 do not have a Gaussian shape and 
we do not show it in this table. 

However, as was discussed in the introduction, the relevant question is whether one can demonstrate with 
confidence that the line is not present in the spectrum of the regions Reg 1-Reg 3 and that the observed 
photon counts are simply statistical fluctuations around a featureless model. Below we argue that this is not 
the case. 

We notice first of all, that the quality of fit (and therefore the conclusion about "chance probability" of 
a given x 2 ) is sensitive to the inclusion or omission of the bins with E > 150 GeV (as comparison of the 
columns (a) and (b) in the Table[T]demonstrates). These bins have extremely low statistics (0 to 2 counts) 
and their large error bars artificially improve the fit quality (as the comparison of the columns for X^iso 
and^g () 20() demonstrates). Inclusion of even larger number of empty bins with E > 200 GeV would further 
improve the quality of fit. At the same time, the change of the interval can change the resulting p-values by 
as much as the order of magnitude (see the columns (c) and (d) in the Table[T|i. 

A way to determine whether a spectral feature is a fluctuation (that is at least partially free from the above- 
described ambiguities) is to add an additional component (Gaussian or other spectrally localized feature) and 
see whether the quality of fit improves. In doing this, we see that for Reg 1, Reg 3 and CENTRAL regions 
the quality of fit improves (by 2 - 3cr). The results are presented in Table [3] 

This method, however, crucially depends on our knowledge of a spectral shape of a feature. As a result, 
we cannot characterize in this way, for example, the region Reg 2, where the feature at 70-100 GeV is 
clearly non-Gaussian. Therefore we also compute the local significance of these features — the Poisson 
probability to observed n counts or more in a given bin, provided that the model predicts A counts. The 
results are presented in the Tableland in Fig. |4] Clearly, such a local significance should be higher than the 
p-value, presented in the table Q] And indeed, in the spectrum of the CENTRAL region we recover an excess 
at 130 GeV at 7.34cr and a 3.6cr excess in the energy bin around 1 10 GeV line. The spectrum of the REG 1 
region demonstrates an excess above the power law fit at energy ~ 115 GeV with the significance 4.63<x 
(Fig.[3a|i. The spectrum of the Reg 2 region (Fig.lJbl formally has a reasonable power law fit in the range 
60-200 GeV (reduced^ 2 = 1.19 for 28 degrees of freedom). However, the distribution of residuals in this 
region shows a broad positive fluctuations in the range 70-120 GeV. Indeed, a power law fit to the subset 
of the data (energies 60-110 GeV) is extremely poor (the reduced^ 2 > 2). However, in the region 110- 
200 GeV the data points have just few counts and therefore large statistical errors reduce the overall^ 2 value 
to acceptable values close to 1 . Finally, in the region REG 3 we observed a line-like feature at 80 GeV (and a 
possible additional feature around 140-150 GeV) with significance around 3cr. Notice that positions of both 
features in Reg 3 coincide with the peaks in the local significance of the Reg 2. We also checked that thus 
determined significance does not depend on the number of bins (the maxima of the lines of different color 
in the panels of Fig. [4] coincide). Notice, that the bin of 20 GeV is much wider than the spectral resolution 
of the Fermi satellite at these energies (about 10 GeV), suggesting that this is not a line-like feature (whose 
width should be 10% of the energy) but rather a power law with a sharp cut-off (as e.g. in 111511 ). 

In order to further clarify the significance of the observed features (and correct for possible "trial fac- 
tors"), we perform 10 4 Monte Carlo realization of the background (power law) model for regions REG 1, 
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Figure 3: Spectra from the Reg 1 region (magenta line), panel |(a)) , Reg 2 (panel [(b)) , REG 3 (panel [(c)) and central region (panel |(d)} , 
The definitions of the regions are shown in Fig. [2] The best fit power law for the central region is shown in black dashed line, the best fit 
diffuse FERMI background model is shown with solid color line. The thickness of the spectra is determined as ±lcr around the central 
value. All spectra together are shown in panel |(e)| The spectrum of Reg 2 with 10 GeV energy binning together with the corresponding 
prediction from the FERMI diffuse background model is given in panel |(f)| f or comparison (see Section l2~3l for details). 



REG 2 and Central (see Fig. [5]). In all panels 99.9% of all realizations are inside the dashed lines. Feature 
at 130 GeV in the CENTRAL region (Fig. |5aJ lies clearly outside these lines, which means that the signifi- 
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(a) REG 1 . A feature at 1 10 GeV is clearly visible 
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(b) REG 2. A broad feature in the region 70 - 100 GeV is 
clearly visible 





(c) REG 3. Possible features are around 80 GeV and around 
150 GeV (compare excess in REG 2 region at these ener- 
gies). Notice that for this region we also added one bin local 
significance. 



(d) CENTRAL. A strong feature around 130 GeV and a pos- 
sible second feature around 110 GeV. 



Figure 4: Local significance of observing fluctuations in 2, 3 and 4 consecutive bins (with 10, 15 and 20 GeV width correspondingly), 
as a function of energy. The significance is defined as the Poisson probability to observe N or more counts, when the model predicts A 
(based on the power law background model). Notice that the significance does not depend on the number of consecutive bins that we 
choose! 



cance of the feature is above 3cr. Feature at 1 10 GeV in Reg 1 (Fig. [5bJ> is also a 3<x deviation^ Notice, that 
in case of the Reg 2 our simulations provide a conservative estimate of significance. The described proce- 
dure does not take into account that in the spectrum of this region several consecutive bins deviate from the 
power law model. This analysis allows us to conclude that, although the spectral feature at the GC is more 
significant (possibly due to the lower predicted background at 130 GeV as compared to lower energies), the 
presence of the spectral features in other regions cannot be ruled out with confidence^ 



7 The discrepancy between the p- value, deduced from Monte Carlo simulations, and p- values, listed in the Tabled is of course not 
surprising — x~ is a global measure of the quality of fit, while the features that we are discussing are localized to several consecutive 
bins. 

8 The regions REG 1, REG 2, REG 3, discussed so far, as well as the regions where the most significant excess is observed in @|, 
are located in the different regions along the Galactic plane. In [2], however, the excess at 130 GeV was claimed to originate from a 
large region, mostly located outside the Galactic plane (its approximate shape is shown in Fig. [5] left panel). Removing the central box 
of 3° X 3° (about 10% of the total field of view) we see that the significance of the feature at 130 GeV drops (green vs. blue data points 
in the right panel in Fig. [9}. Therefore, we conclude that the feature at 130 GeV is related to the region close to the Galactic center, 
without inclusion of this region the feature becomes insignificant (a similar conclusion was reached in (3),0I). 
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(b) Region REG 1. Result of 10 4 Monte 
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ture at 1 10 GeV is more than 3<x. 
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each deviation is above 2.3<x. 



Figure 5: Monte Carlo realizations of the power law models in the region 50-200 GeV for CENTRAL (left), REG 1 (center) and REG 2 
(right) regions. Number of realization is 10 4 . Shaded regions show minimum and maximum number of counts in all the Monte Carlo 
realizations. 99.9% of all realizations lie between dashed lines (in panel |(c)| dashed-dotted line show 99% of all models). Data points 
outside dashed lines are at least 3<x deviations from the power law model. Two data points in the panel |(c)| deviate from the power law 
by 2.3<r. 

2.3. Test statistics maps based on Fermi diffuse background model 

One possible interpretation of the above results is that by exploring regions outside the GC one increases 
the "trial factors" by 10° x 60° j(n x 3° x 3°) ~ 21 (and therefore p-vales of the features, obtained above, 
should be further penalized by this factor). This would be the case if spectra in all 3° circles around the 
Galactic plane were completely arbitrary, not correlated with each other. This is however, not true in our 
case, as the diffuse emission in the Galactic Plane is described by the smooth model. 

To study both spectral and spatial variations of emission around the Galactic plane, we model 50° X 50° 
region (in J2000 coordinates) around the GC by three components — standard galactic(centered at GC) and 
extragalactic Fermi backgrounds^ and a model of a "line", that we try to add in each pixel in order to improve 
the fit (pixel size is 0.5° x 0.5°). The "line" model has Gaussian spatial profile with 3° FWHM and Gaussian 
spectral profile with 5 GeV dispersion, centered correspondingly at 80 GeV, 110 GeV, 130 GeV"F"l We fit 
the data with the diffuse models and then consider the change of the log-likelihoods of the fit when adding 
the Gaussian component at these three energies. We build the map of TS values (the difference of log- 
likelihoods of two models) using gttsmap tool. The normalizations of the Fermi diffuse backgrounds, as 
well the normalization of the region model were allowed to vary during the fit procedure. The test-statistics 
values of the 3°-region model fits are shown in Figs. [6] The significance of the model can be estimated as 

VTScr. 

The test-statistical value of the line centered at CENTRAL region is about 21, that corresponds to signifi- 
cance ~ 4.5<r. Beside the CENTRAL region Figs.|6]exhibit additional regions with the excesses at 70-90 GeV 
(Fig.[6cji and 100-120 GeV (Fig.|6b]i. Notice that, contrary to the results of 0], the test statistics of the CEN- 
TRAL region does not improve significantly when adding a Gaussian at 110 GeV, i.e. photon distribution 
in this energy bin strongly deviates from 3° FWHM spatial Gaussian profile. The highest test statistics 
TS = 32.3 is found in the region Reg 2 at 80 GeV (see Table[3]for details). 



Given by gal_2yearp7v6_v0 . f its and iso_p7v6clean.txt templates, see e.g. 

|http : lit ermi . gsf c ■ nasa ■ gov/ ssc/data/ access/la t/BackgroundMode fs". html| 

10 We have repeated the same exercise with the Gaussian line, having 2 GeV FWHM and obtained similar results. 
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(a) Energy range 120-140 GeV. The maximal TS s21.2 is located at the pixel with (/, b) a (-1.5°, -0.75°). 




.2 4 6 8 10 12 14- 16- 18 

(b) Energy range 100-120 GeV. The maximal TS in REG 1 is TSx 15.0 is located at the pixel with (/, b) * 
(18.1,0.2) and in the region REG 2 the maximal TS is a 17.3, located at the position (l,b) = (337.1,-0.2) 
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(c) Energy range 70-90 GeV. The maximal TS »32.3 is located at the pixel with (/, b) as (336.0, -0.2). 

Figure 6: Test-statistics (TS) maps when in addition to Galactic and extragalactic Fermi backgrounds we add a model of the 3°-region 
centered at this point. The latter was taken to has Gaussian spatial profile with 3° FWHM and Gaussian spectral profile with 5 GeV 
dispersion, centered correspondingly at 80 GeV, 110 GeV, 130 GeV. 
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Region (Z°, b°) 


Energy of the line 




Reg 1 (18.72°, 0.57°) 


HOGeV 


14.99 


Reg 2 (335.76°, 1.255°) 


80GeV 
HOGeV 


32.3 
17.26 


Reg 3 (12.22° ,0.20°) 


80GeV 


2.74 


Central (359°,-0.7°) 


130 GeV 
HOGeV 


21.22 
4.5 



Table 4: Change of the test statistics (TS) when fitting spectral and spatial distribution of photons to a smooth background model and 
to the same model with an additional line (see Section ff) The line has Gaussian spatial profile with 3° FWHM and Gaussian spectral 
profile with 5 GeV dispersion on top of Fermi diffuse background models. The TS maps are shown in Fig. [6] The TS max column shows 
the maximum difference of likelihoods over all pixels within a given region. 

3. Discussion and conclusions 

The first version of this paper (1205.4700vll demonstrated that spectra of several regions (e.g. Reg 1, 
REG 2, Reg 3) are not featureless at energies E > 50 GeV, but contain excesses similar in significance 
to the spectral feature in the GC region. In particular, version 1 found excess in 110 GeV bin from the 
REG 1 (a pair of lines at ~ 110 and ~ 130 GeV from the GC region was later discussed in |]4|]). Our 
present analysis recovers the previous results with high significance and our main conclusion remains intact 
- with the current data it is impossible to rule out the possibility that the observed spectral features are 
actual signals (of astrophysical or instrumental origin), rather than statistical fluctuations of the smooth 
backgrounds. Clearly, until such a possibility is ruled out, the DM interpretation of the emission from the 
GC remains dubious. 

Refs. ll4 ll 811 also analyzed spatial distribution of residuals and significance of different regions around 
the Galactic plane. The analysis, presented here, and the interpretation of the significance of the spectral 
features is different due to the different statistical approaches adopted. 

One way to detect spacial "hot spots" with significant spectral features is just to look for the regions 
where the best fit to a smooth (power law) model is unacceptable (i.e. the reduced^ 2 » 1). This method 
was discussed e.g. in [20, 19]. The results of such fits for our regions were presented in TableQ] This method 
however, is prone to a number of uncertainties, related to the choice of fitting interval and to the fact that y 2 
is a global measure of fit, where the significance of local spectral features can be "blurred" by many bins 
with large errors bars. Indeed, as demonstrated in Section l2~2l the local significance of the features, found in 
regions Reg 1, Reg 2, and Reg 3 is higher than predicted by this method (3cr and above). 

A common way to determine the significance of these hot spots is to multiply the resulting p-values of 
the spectral features by the number of regions that were used in the search for the "hot spots" (so called "trial 
factors", see e.g. ||4[]). The number of regions is, however, totally subjective and the resulting significance 
can be artificially lowered (or increased). We believe that a better way to account for spatial distribution 
of residuals is to find a smooth background model that fits the data in (most of the) spatial regions and 
deviation from which are normally distributed (modulo some localized spots) (a similar approach was used 
in J3l). We used the galactic diffuse background by the Fermi collaboration as such a spatial model. We see 
(Figs. [6j that this model correctly fits the data "almost everywhere" apart from several localized spots where 
it underpredicts the number of photons and where additional Gaussian strongly improves the quality of fit. 
The histogram (Fig. |7]l clearly shows that the distribution of pixels with TS < 8 is consistent with the \ 2 
distribution with 1 degree of freedom (as it should be for statistical fluctuations). However, there are also 
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few high TS outliers, corresponding to the hot spots in Figs. l6aH6ctF^l These regions are the GC and several 
other regions (Reg 1 and Reg 2) that we have previously identified in Figs.[TJ|2] 

Looking at Figs. [3] one notices that in these regions Fermi model systematically (in many consecutive 
bins) underpredicts the data (this is also true for the GC region){3 As Fermi diffuse model is determined 
globally, one cannot improve the quality of fit in these spots without worsening the overall fit. In all these 
regions the signal is consistent with the smooth background model plus additional sources of emission with 
sharp energy spectra localized in the Galactic plane. 

Indeed, let us compare two most prominent regions (CENTRAL and Reg 2) (Figs. [3b1 and [3dTi . We can 
interpret their spectra as a combination of spatially smooth (Fermi) background model plus a localized excess 
modeled in each region by two components: a low-energy (20-80 GeV) power law and a feature at energies 
around - 70-80 GeV for Reg 2 and 1 10- 130 GeV for CENTRAL regions The count rates observed in these 
features (number of counts in excess of the background model) are approximately the same, the significance 
of the feature is thus determined only by the predicted background value. This value of the predicted signal 
decreases with energy, therefore the most energetic (130 GeV) feature, having the count rate close to that of 
features at other energies, would have a higher significance against very low background at that high energy. 

An additional evidence in favor of the "local source" hypothesis is provided by the histogram of the 
number of point sources from the 2 year Fermi catalog along the Galactic Plane (shown in Fig. [8}. One sees 
that the regions CENTRAL, Reg 2 and Reg 1 correspond to the local maxima in the concentration of point 
sources in the Galactic plane. This picture is consistent with having some additional localized sources of 
emission in some regions along the Galactic plane. 

In conclusion, at the moment there are two possible interpretations of the observed features. In one of 
them the 130 GeV feature is assumed to be a unique statistically significant feature (probably explained by 
DM annihilation/decay) see e.g. JUS], while all other features are treated as fluctuations of the powerlaw 
background. The other possible interpretation is that all the features at 80 GeV, 110 GeV, 130 GeV observed 
in the spectra of different parts of the Galactic plane have a similar origin. Our analysis shows that the present 
data is consistent with "many features - many sources" interpretation. If confirmed, such interpretation 
would be hard to explain in terms of Dark Matter annihilation (or decay), it would probably require an 
astrophysical explanation, like, for example the one of lfl5ll . in which narrow spectral features (located at 
different energies for different systems) are produced by emission from ultra-relativistic pulsar wind. 

If this interpretation will be confirmed by the future data, this means that a line-like feature in gamma 
ray spectrum should not necessarily be "a smoking gun" of dark matter annihilation (or decay), as it is often 
assumed. In Refs. J2, 0. fl8ll a procedure of data analysis is described in which only the most significant 
(located in the GC) feature survives while all other regions have their significances below 3<x. However, 
to speak seriously about the DM origin of a signal, one has to show that no other interpretation is possible 
(for example, that the other features may not have physical origin). We believe that for the moment the 
data does not allow to distinguish reliably between these two interpretations. Additional observations with 
with HESS-2, Gamma-400 and CTA will probably be required in order to check these models (see [22]), but 
until that the DM interpretation of 130 GeV feature and treating all other features in a different way looks 
dubious. 
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the dashed line and as a result there are not pixels with TS > 8 (for the number of pixels that we have in our regions). However, if 
high-TS outliers for CENTRAL, REG 2, REG 1 , etc. are not caused by statistical fluctuations only, one should see a tail of large TS (as 
indeed seen in this Figure). 
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Figure 8: Number of 2FGL sources at different positions along the Galactic plane (\b\ < 3°). Positions of regions CENTRAL, REG 2 
and REG 1 are clearly visible as local maxima in the number of sources. 
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Figure 9: Left panel: region that approximately coincides with the "Region 3" of | 2] and reproduces its spectrum. Right panel: 
Spectrum of this region (blue) with the feature at 130 GeV clearly visible. The same spectrum extracted from the region with removed 
central 3° X 3° box (green) and decreased significance of 130 GeV feature. 
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